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We study the low-energy properties of a Hubbard chain of finite size Nc connected to two 
noninteracting leads using the numerical renormalization group (NRG) method. The results 
obtained for Nc ~ 3 and 4 show that the low-lying eigenstates have one-to-one correspondence 
with the free quasi-particle excitations of a local Fermi liquid. It enables us to determine the 
transport coefficients from the fixed-point Hamiltonian. At half-filling, the conductance for even 
Nc decreases exponentially with increasing U showing a tendency towards the development 
of a Mott-Hubbard gap. In contrast, for odd Nc, the Fermi-liquid nature of the low-energy 
states assures perfect transmission through the Kondo resonance. Our formulation to deduce 
the conductance from the fixed-point energy levels can be applied to various types of interacting 
systems. 
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' 1. Introduction 

Electron transport through finite systems, such as 
, quantum dots, quantum wires, and atomic chains of 
' nanoscale, is a subject of much current interest. In these 
, systems, a number of phenomena have been predicted 

■ theoretically and some have already been successfully ob- 
, served. The Kondo effect in quantum dots is one such 

■ example. Furthermore, recent experimental develop- 
ments make it possible to examine the interplay of var- 
ious effects which have previously only been studied in 

, different fields of physics. For instance, in quantum dots, 

■ the interplay of the Aharanov-Bohm, Fano, Josephson, 
, and Kondo effects under equilibrium and nonequilib- 

■ rium situations have been studying intensively.^"^ The 
, Luttinger-liquid behavior in quantum wires has also been 

■ an active field of investigation,®' ^ and numerical develop- 
ments have been reported recently for spinless fermions 
on a lattice. 

We have previously considered the transport proper- 

■ ties of a finite Hubbard chain of size Nc (=1, 2, 3, ... ) 
I connected to two noninteracting leads, which is illus- 

■ trated in Fig. 1 , as a model for a series of quantum dots 
] and materials on a nanometer scale. We have discussed 
I an even-odd oscillatory behavior of the dc conductance 

at half-filling based on a perturbation expansion in pow- 
ers of the Coulomb repulsion L/.^^"^^ For even Nc, the 
conductance decreases with increasing U showing a ten- 
dency towards the development of a Mott-Hubbard insu- 
lator gap.^^"^'' The conductance deduced from the order 
self-energy was qualitatively correct. However, a more 
accurate treatment is needed to deal quantitatively with 
the large U regime, which is one of the aims of this paper. 

In contrast, for odd Nc, the transmission probabil- 
ity through the Hubbard chain reaches the unitary-limit 
value at r = when the system has the electron-hole 
and inversion symmetries. Physically there is a Kondo 
resonance situated at the Fermi level which enables the 



perfect transmission take place for any value of U. The 
proof was given by taking all contributions in powers of 
U formally into account for the Nc x Nc matrix self- 
energy S(ti^). The assumption we have made is; ReS(O) 
is not singular and Im5](0) = at T = 0.^'^ The local 
Fermi-liquid state satisfies this assumption, and his- 
torically the same assumption has been made by Langer 
and Ambegaokar in the derivation of the Friedel sum 
rule for interacting electrons. In the case of the sin- 
gle Anderson impurity corresponding to Nc = 1, the 
perturbation theory in powers of U describes the low- 
energy Fermi-liquid behavior, and is consistent with 
the exact Bethe ansatz^^^^o ^^^^ NRG^^^^^ results. The 
perturbation expansion in U works because the contri- 
butions from the low-energy processes, in which elec- 
trons hop into the reservoirs and away from the impu- 
rity, are included in the noninteracting Green's functions 
via the hybridization energy scale F^ and T/j, where 
^L/R = ''^P'^'l/r- The perturbation expansion is also con- 
vergent for Nc > l}'^^^'^ It would be of interest to have 
these results confirmed by a non-perturbative technique, 
which is the second aim of the work presented in this 
paper. 

To tackle both these problems we apply the non- 
perturbative NRG approach to the low energy physics of 
the Hubbard chain, connected to non-interacting leads. 
In doing so we go beyond the earlier low order perturba- 
tional results for the even site chains, and we also derive 
a Fermi liquid picture for the odd site chains, without 
making the assumptions implicit in the perturbation the- 
ory. The NRG method has been applied successfully to 
the quantum dots for Nc = 1 and 2.^'*'^''^ In the present 
work we have applied the NRG method for the Hubbard 
chain with size Nc = 3 and 4, which seem to capture the 
essence of the even- and odd-size chains, respectively. 
The results show that the low-lying energy states have a 
one-to-one correspondence with the quasi-particles exci- 
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Fig. 1. Schematic picture of system. 



tations of the local Fermi liquid. It assures the validity of 
the Fermi-liquid description at low energies. Specifically, 
for odd Nc, a number of the noninteracting sites are 
needed to be taken into account to reach the fixed point 
that describes the physics below the Kondo temperature 
Tk- The fixed point has much information about the 
low-temperature properties, and one can deduce the pa- 
rameters such as Tk and Wilson ratio R from the flow 
of the eigenvalues.^^' In the present paper, we provide 
a formulation to determine the conductance for even Nq 
at T = 0. For large U, the NRG method improves the 
perturbation results, and the conductance determined 
from the fixed-point energy levels decreases exponentially 
with increasing U . 

In §2, we deduce expressions of the ground-state prop- 
erties in terms of the Green's function. In §3, we de- 
scribe the formulation to deduce the conductance from 
the fixed-point Hamiltonian. In §4, we show the NRG 
results. In §5, discussion and summary are given. 

2. Model and Formulation 

Wc consider a Hubbard chain of a finite size Nc 
situated at the center, which is attached to two non- 
interacting leads at the left {L) and the right [R). The 
complete Hamiltonian is given by 



W = Hd + Hu + Wmix + Wlead , 



(1) 



Nc-1 

E E {4A+la + 4+lAa) 
i=l a 



Nc 

+ Ei«<i 

i=l 



U 

]ndi, 



(2) 



(3) 



i=l 



+ ^rY.{ '^R^'^Nca + 4c,>fl, ) ' (4) 

a 

ij—L,Ii k(j 

where dia annihilates an electron with spin a at site i, 

and Hdi = J2a4.a^irT- ^^^^ ^' -^)' 

operator cj^^^^ creates an electron with energy e^i, corre- 
sponding to an one-particle state (j)kv{r). The hopping 
matrix elements Vj^ and Vj^ connect the chain and leads. 
At the interfaces, a linear combination of the conduction 



electrons 



^kvirv) mixed with the elec- 



trons at i = 1 or Nc (as illustrated in Fig. 1), where 
denotes the position at the interface in the lead side. 



For this system, the Green's function is defined by 



(6) 



where /3 = 1/T, dj^ir) = e^^d 



-tH 



, and (• • • ) de- 
notes the thermal average Tr [e~^^ • • •] /Tre~^^. We 
use units h = 1. The corresponding retarded function, 
G^j,{lj) = Gjji{uj + iO+), is obtained via the analytic 
continuation. Since the interaction U is finite only for 
the electrons in the chain at the center, the Dyson equa- 
tion is written in the form 

Nc 

Giiiz) = G^iz) + G°Az)i:i,Az)GMz) . (7) 

i'j'=l 

Here G'^j{z) is the unperturbed Green's function corre- 
sponding to Ho = H-d + Hiead + '^mix) and Sjj(2:) is the 

self-energy correction due to Hu- Note that Gij{z) = 

Gji{z), because of the time reversal symmetry of H.. The 
Dyson equation can be rewritten using Nc x Nc matrices 
G{z) = {Gij{z)} and T,{z) = {T'ijiz)} as 



{G{z)y 



where 



zl 



~t 



Vxnix(^) - S(^) , (8) 



^mix('2^) — 















-t 



(9) 











'"rSr{z). 



(10) 



and g+{uj) = Efe |<^,.ct(^i-)| /(t^ - ffe,/ + «0+) is the 
Green's function at interface of the isolated lead. In the 
present study, we assume that the density of states is a 
constant, gt{^) = — iTT/Cy, for small u>. Then the energy 
scale of the level-broadening becomes = np,yvl, and 
it determines the two non-zero elements of V^;^(a;). 

2.1 Ground-state properties 

If the ground state has a property ImS"'"(0) = at 
T = 0, the damping of the excitations at the Fermi level 
vanishes. Then the effective Hamiltonian defined by 

H^^ = iJ^-FReS+(0) 



(11) 

plays a central role on the ground-state properties. It 
determines the renormalized hopping matrix elements 
H'^ = {—tij }, and also the value of the Green's func- 
tion at the Fermi level {G+(0)}"^ = K{0) - V+-^{0), 
where 

K{uj) = Ljl- . (12) 

The determinant of the matrix {G~''(0)} ^ is related to 
the scattering matrix, and can be rewritten in the fol- 
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lowing form by expanding the first and A^c-th columns, 
det{C?+(0)}"^ = \-TLTiidetK^f^''{0) + detK{0) 



Tr detKjvciVcW 
(13) 



Here (0) is a {Nc - 1) x {Nc - 1) derived from K{0) 
by deleting the i-th row and the j-th column. Similarly, 
Kff ^^(0) is a {Nc-2) x {Nc-2) matrix obtained from 
K{0) by deleting the first and the Nc-th rows, and the 
first and the A^c-th columns. At T = 0, the dc conduc- 
tance is determined by the Green's function which con- 
nects the two leads, gj^^ = {2e'^ / h) AT rT l |G;;^c,i(0)|^. 
It can also be expressed in terms of the scattering matrix 



3nc 



h 



4rr,r 



LJ- R 



[detKiNciO]] 



det{G+(0)}' 



(14) 



Furthermore at T = 0, the charge displacement can be 
determined by the Priedel sum rule, 



- Im log det{G+(0)}" 



(15) 



Particularly for the constant density of states, eq. 
(15) corresponds to the charge displacement defined by 
ANioi = Recently, a related formulation 

which takes into account the self-energy corrections us- 
ing the effective Hamiltonian has also been applied to a 
finite ring with a magnetic flux.^® 

2.2 Conductance at half- filling 

Specifically at half-filling ej, = —U/2, the matrix ele- 
ments of K{0) = {tij } become zero for i and j belonging 

to the same sublattice, i.e., \i — j\ = 0, 2, 4, Thus in 

this case, K{0) has a checkered structure, and it causes 
the even-odd dependence on the number of the interact- 
ing sites Nc- 

For even Nc (= 2M), eq. (14) can be rewritten in the 
form 



_ 2e^ 

92M — u 



~2 _ _ det JC(0) 



(16) 



(17) 



where we have used the relations which can 
be deduced from the checkered structure of 
K{Q): detKMcNci'^) = 0, detKii(O) = 0, and 
[detKAr,,i(0)]^ = -deiK{0) deiK^^^'^{Q). Specifi- 
cally, for free electrons at U — 0, K{{)) is given simply 
by -H%, and eq. (17) yields v^ = t^. 

For odd Nc (= 2M-|- 1), the dc conductance, eq. (14), 
can be expressed as 



92M+\ 



2e' 



R 



h [{TL+VR)l2f 
where f z, = AFz,, f ^ = F/j / A, and 



A 



/ detail (0) 
deiKNcNci^) 



(18) 



(19) 
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Fig. 2. Schematic pictures of discretized Hamiltonian for the 
NRG approach for (a) iVc = 3 and (b) AT^ = 4 



Here we have used the properties, detii'(O) = 

0, detKfi^^^(O) ^ 0, and [det KAr^,i(0) ]^ = 
detifii(O) detlf jVcJVc(O)- Especially, if the system has 
the inversion symmetry = F/j in addition to the 
electron-hole symmetry, the parameter defined in eq. (19) 
becomes A = 1, and then the perfect transmission occurs, 
52M+1 — 2e^//i, for any value of M and U. 

3. Fixed-point Hamiltonian and conductance 

In the NRG method, the conduction band can be mod- 
elled by a linear chain as shown in Fig. 2, via a stan- 
dard procedure of logarithmic discretization.^^' Then, 
to capture the low-energy behavior correctly, we use a 
sequence of the Hamiltonian defined by 

= k^^-'y^^Ha + nv + H^^ + H^^y), (20) 



+ Vr ^ (/o.fl.C?JVc,^ 



) , (21) 



(N) 
lead 



D 



1 + 1/ A V- v-v- 



N-l 
v=L,R <7 n=0 



^ ( fn+lyua- fn,va fn,va fn+l,i/a ) ' (^2) 

where D is the half-width of the conduction band. The 
hopping matrix elements and ^„ are defined by 



2-DF^Aa 



1 1 + 1/A 

2 1 - 1/A 



1- l/A"+i 



v/l-l/A2»+Vl-l/A2^ 



n+3 



log A, (23) 



(24) 



The factor is needed to compare the discretized model 
with the original Hamiltonian eq. (1) precisely, and it 
behaves as Aa ^ 1 in the continuum limit A ]^ 22,30 
The low-lying energy states of the original Hamiltonian 
H can be deduced from those of A~^^~^^/'^Hn for large 
N. 
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In the following we concentrate on the case = F/j 
(= r), where the couplings to the two leads are symmet- 
ric Vj^ = Vj^ (= v). In the discretized Hamiltonian Hn in 
eq. (20), the matrix elements t and v has multiplied by 
7Y(^-i)/2 j\g shown in the next section, for large A'' the 
low-lying energy states of the many-body Hamiltonian 
Hj\f converge to the states which have; one-to-one corre- 
spondence to the quasi-particles of a local Fermi liquid. It 
enables us to deduce the matrix elements of H^, which 
are defined in Eq. (11), from the NRG spectrum. 

At the fixed point, the low-energy spectrum of the 
many-body Hamiltonian can be reproduced by the 
one-particle Hamiltonian consisting of and the two 
finite leads; 

<^ = A(^-)/^ ( Hf + H^,. + H W ) , (25) 
~ S«=i '^L^i-- describes the free 



where 



,JVc 

quasi-particles in the cluster with Nc + 2(7V -I- 1) sites, 
and the corresponding Green's function can be written 
as 



{G,p(a;)}-^ ^ 

^(JV-l)/2 LA-(iV-l)/2 ;^ 



rj-eff 



(26) 



Here we have not included the renormalization factor 

d'S/duj, because at T = it docs not affect the dc con- 
ductance and charge displacement defined in eqs. (14) 
and (15). An eigenvalue s* of Hqp^ satisfies the equation 
det{Gqp(e*)}~' = 0, which can be written in a form 
similar to eq. (13), 

- [detKuM +det KnoNcM] [v''A^''-^^/^gj,{e*)] 

= , (27) 

where lon = e* A^(^^^)/^. The Green's fmiction gjv(a^) 
is introduced for an isolated lead with -|- 1 sites, and 
is defined with respect to the interface n = 0. It can be 
expressed as g^{uj) = YlZ=a l'^™(0)P/(^ ~ ^rn), where 
€m and <^m('^) are the eigenvalue and eigenstate for the 
isolated lead. In the electron-hole symmetric case, Eq. 
(27) can be simplified by using the properties described 
in Sec. 2.2, as follows. 

For even Nc (= 2M) and in the Umit of large A'', Eq. 
(27) yields 

2 detK{0) 



hm A(^- 

N—>oo 



l)/2 



detK^f^^iO) 



(28) 

Thus, the parameter Vq defined in eq. (17) can be related 
to the eigenvalue e*, as 



r2 



V 

TD 



N 



lim DA(^-i)/2gjv(£*) 



(29) 



The prefactor in the right-hand side can also be written 

as v'^/{TD) = lAt^j-K by using eq. (23). Note the func- 
tion, limjv-»oo A.(^~^^/2gjy(^(^^^ depends on whether A'' is 
even or odd. For even Nc, the fixed-point eigenvalue e* 



depends on U jt and F/t, and the fixed-point Hamilto- 
nian can be written in the form 



(30) 



where A^qp = Nc/2 + N ■ 



1. One can determine w^/F^ 
by substituting the value of deduced from the NRG 
results into eq. (29), and then the dc conductance can 
be obtained from eq. (16). Our formulation to deduce 
the conductance is analogous to the method used for the 
asymmetric Anderson impurity to determine the local 
charge from the fixed-point eigenvalue and Friedel sum 
rule.^^ 

The situation is quite different for odd Nc (= 2M-I-1). 
In this case, eq. (27) yields two separate branches of low- 
energy states for large A'', 



and 



lim A(^-i)/2g^(£*) = , 

AT— »(X) 



lim A(^-i)/2g^(e*) = 



(31) 



(32) 



These can be deduced from the behavior of the de- 
terminants for small wjv, det X^'^^'^ (wjv) oc lon and 
det K{u)n) oc u)n- These two branches imply that the 
eigenvalue e* does not depend on [/, t and F. Specifi- 
cally, eq. (31) corresponds to an isolated lead consisting 
of A'' sites starting from n = 1 and ending at n = AT, 
while eq. (32) corresponds to another lead with size A^+1 
that includes the site n = 0. One way to interpret the 
fixed point is along the lines of the original work of Wil- 
son, as a strong coupling fixed point, such that a single 
site at n = is removed from one of the leads to join 
the interacting sites and to form a singlet ground state 
for this cluster consisting of A'c' + 1 sites. However, here 
we derived eqs. (31) and (32) for the connected chain of 
Nc + 2{N +1) sites and have interpreted the fixed-point 
using the Hamiltonian h'^^ , which is defined in eq. (25) 
and can be diagonalized as 



Oct "0(7 



(33) 

where el = Q and AT^p = {Nc - l)/2 + N + 1. This 
Hamiltonian links directly to the Fermi-liquid behavior, 
as the quasi-particles defined in this connected chain are 
in one-to-one correspondence with the single-particle ex- 
citations of the non- interacting system {U = 0). This 
leads to more natural description for the ground state 
properties than the strong coupling interpretation which 
involves breaking the chain by effectively removing two 
sites. As we show in the next section, the low-lying 
energy states of the many-body Hamiltonian Hm repro- 
duce the energy spectrum determined by eqs. (31) and 
(32). The same behavior was seen in the single impu- 
rity case.^^'^^ Our numerical results confirm the Fermi- 
liquid behavior for Nc = 3, and justify that the assump- 
tions made in deducing the unitary-limit transport re- 
sult, g2M+i — 2e^//i, for the case A^c = 3. 
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Fig. 3. Low-lying energy levels of Hf^/D for Nc = 4 as a function of odd A'' (up to 29) for several values U/{2-Kt): (a) 0.0, (B) 0.5, (c) 
1.0, and (d) 1.5. Here = — 1//2, V /t = 0.12, t/D = 0.1 and A = 6.0. The eigenvalues are measured from the ground-state energy for 
each N . The label {21 , 2S) corresponds to the total axial charge I and spin S. For A > 7, the levels approach to the fixed point values. 



4. NRG results 

In the electron-hole symmetric case, the Hamiltonian 
Hn has a global SU{2) symmetry of the axial charge, '^^ 
which is specified by the generators 

Nc 
i=l 



N ^ 

V V - f f t f + f - 1 



u=L,Rn=0 



(34) 



AT 



u=L,B.n=0 



N 



(35) 



J_ = ^(-l)Hxd.T + E E(-l)'"^/n,^i/».^T 



iy=L,Rn=0 



(36) 

Here 6',i.L = —n for the left lead, and On,R. = A'^c + ft + 1 
for the right lead, so that the factor (— l)^"-' becomes 
+ 1 or —1 depending on whether the site labeled by 
(n, v) is in an even or odd sublattice. The z compo- 
nent of the axial charge corresponds to the total charge, 
Q — 2Iz- Furthermore, the operators and I± satisfy 



the commutation relations identical to those of the total 
spin operators Sz and S±. Thus, using the symmetry of 
S't7(2)spin X S'[/(2)axiai, the eigenstates can be classified 
according to the quantum numbers for the operators Sz, 
S^, %, and P = ll + {T+T^ + T-T+)/2, as 



HN\I,Iz,S,Sz\r)N ^ ENjs,r\I,Iz,S,Sz;r)N ■ (37) 

The use of the SU{2)spi-a x S'C/(2)axiai symmetry has 
a great numerical advantage [see also, Appendix]. One 
can save the eigenstates to be retained in the process 
of the NRG iteration. Particularly in the multi-channel 
systems, such as the one we are considering (2-channel 
in our case), the number of low-energy states to be re- 
tained increases exponentially with the number of the 
channels. Thus the reduction of the Hilbert space helps 
to improve the numerical accuracy. After using the two 
SU{2) symmetries, we have retained typically 1000 low- 
energy states in the NRG calculations for the Hubbard 
model with Nc = 3 and 4. Under this condition the 
truncation occurs first at = 2, namely Hn has been 
diagonalized exactly up to A = 1 where the total num- 
ber of the sites in the cluster is 7 and 8 for Ac = 3 and 
4, respectively, as shown in Fig. 2. Since two new sites 
at n = A + 1 are included in each step of the recursive 
procedure, the dimension of Hilbert space becomes typ- 
ically 4^ larger than the number of the states retained 
in the previous step apart from some reductions due to 
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Fig. 4. The {/-dependence of e* and fo'' the 4-site Hubbard 
model for odd N. For large U, these two approach to a single 
value 0.5312, which corresponds to the lowest pole of the Green's 
function of the isolated lead, linijv^cx) ^^^~^^^^SNi'^)' °'1<^ 
N. 



the symmetries. To overcome the influence of the trun- 
cation, which starts not so far from the interacting re- 
gion, we have concentrated on the case in which F/i is 
small (~ 0.1) and have used a rather large value for the 
discretization parameter A = 6.0. The hopping matrix 
between the interacting sites has taken to he t/D = 0.1 
in the calculations. 

The use of an additional inversion symmetry T^^ = 
does the change the total number of the basis states 
to be retained in the interacting case U 0, although 
each subspace labeled by / and S can be divided up 
into two segments. The inversion symmetry can be em- 
ployed by introducing the bonding and antibonding or- 
bits a, 



n,it,cr 



{fn,R,a^fnX,a)l Howcvcr, for even Nc 
these orbits make the axial charge nonlocal, while the lo- 
cality is preserved for odd Nc- We have also performed 
the calculations using these orbits for odd Nq- It makes 
the computer time somewhat shorter, but is not essential 
for improving the numerical accuracy. 

4.1 Results for Nc ^ 4, 

We first of all consider the Hubbard chain of the size 
Nc = 4, which can be regarded as a simplest case for 
an even interacting chain. In Fig. 3, the low-lying energy 
levels of i?jv for V /t = 0.12 are plotted as a function of 
odd N, where iV -I- 1 is even, for several values of the 
Coulomb interaction: U/{2nt) — 0.0, 0.5, 1.0, and 1.5. 
The eigenvalues are measured from the ground-state en- 



El 



2 St 
el 4 

3 el 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 1.6 1.8 



U/(27Ct) 



Table I. Low-lying fixed-point eigenvalues of Hj^ for A^hub = 4, 
for odd N. 



Fig. 5. The [/-dependence of (a) and (b) the dc conductance 
jjY^ for the 4-site Hubbard model for several values of the hy- 
bridization: r/t = 0.04, 0.07973, and 0.12. The parameter is 
deduced form the fixed-point eigenvalue e^ by using eq. (29). 



ergy for each N. The flow of the energy levels converges 
to the fixed-point values for > 10 in these examples. 
In the case of even Nc, the fixed-point eigenvalues de- 
pend on U and T. The numerical results for the many- 
body low-lying states can be compared with the quasi- 
particle eigenstates defined in eq. (30). We found that 
the many-body eigenvalues E* of Hn shown in Fig. 3 
have correspondence to the excited states described by 
the quasi-particles with the energy e^*, as summarized 
in Table I. The first two excitation energies of Hn, i.e., 
El and i?2, determine the two one-particle energies e* 
and respectively. The many-body eigenvalues above 
these two, E^, E^, . . ., agree well with those calculated 
from el and £2 with the assignments given in the table 
I. We have confirmed that the two quantum numbers, 
the total axial charge / and spin 5, for the low-energy 
eigenstates are consistent with these assignments. This 
feature of the low-lying energy states is similar to that 
for the Kondo and Anderson models, and means that 
the low-energy properties can be described by the local 
Fermi-liquid theory. Specifically, for even Nc, the fixed- 
point Hamiltonian Hqp ^ defined in eq. (25) can be sep- 
arated into a couple of the chains when the system has 
an inversion symmetry, and each of the chains can be 
mapped onto a noninteracting version of the asymmet- 
ric Anderson model. ^'^ In Fig. 4, the U dependence of 
el and £3 ^tre shown for the parameter set; T/t = 0.12, 
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Fig. 6. Low-lying energy levels of Hf^ /D for Nc = 3 as a function of odd A'' (up to 59) for several values U/{2nt): (a) 0.0, (b) 0.5, (c) 
1.0, and (d) 1.5. Here e^j = —(7/2, F/t = 0.12, t/D = 0.1 and A = 6.0. The eigenvalues are measured from the ground-state energy for 
each N. The label {21, 2S) corresponds to the total axial charge I and spin 5. The size of the lead, TV*, required to get the energy 
levels closed to the fixed-point values increases with U. 



t/D ~ 0.1, and A = 6. The difference between £* and 
e\ decreases with increasing U . This tendency hnks with 
the behavior of the other excitation energies shown in 
Fig. 3, namely i?3, £'4, and £'5 become close to each 
other with increasing U . For large C/, both £\ and £3 ap- 
proach to 0.5312 which corresponds to the smallest pole 
of limAT^^oo A(^~^^/^g'^((jj) for odd N. It means that in 
the limit of C/ — > 00 the low-energy states are determined 
by those of the isolated leads. 

Substituting the value of e\ or into eq. (29) and tak- 
ing A = 99 for the Green's function of the isolated lead, 
we obtain the parameter vc defined in eq. (17), and then 
determine the dc conductance via eq. (16). The results 
are plotted against U / {2'Kt) in Fig. 5 for several values of 
hybridization; T/t = 0.04, 0.07973, and 0.12. The value of 
Vc determined from e\ agree with that determined from 
£2- Furthermore, we have also calculated vc using the 
fixed-point eigenvalues for even A and the Green's func- 
tion A(^~^)/^g^(a;) for A = 100. The result again agrees 
with that deduced from the data for odd A. The param- 
eter Vc increases with the Coulomb interaction U reflect- 
ing the behavior of the fixed-point eigenvalue shown in 
Fig. 4. The results for the conductance are plotted in 
Fig. 5 (b). When the coupling with leads F increases, the 
conductance also increases in the parameter region we 
have examined. For large U , the conductance decreases 



exponentially with increasing U. This can be understood 
as a tendency towards the development of a Mott insu- 
lator gap. Namely, for large even Ac, the conductance 
is expected to show a behavior gj^^ oc e"^'^/'', where 
^ ~ ?iUi?/Agap is a correlation length determined by the 
Hubbard gap Agap and Fermi velocity vp- This is because 
Agap cx U for large U , as can be seen in the Bethe ansatz 
results for one-dimensional Hubbard model. Our pre- 
vious results obtained with the 2nd order perturbation 
theory in U are valid qualitatively."'^^ However, the 2nd- 
order perturbation theory fails to reproduce the correct 
exponential dependence for large C/, and it has now been 
corrected with the non-perturbative NRG technique. 

4.2 Results for Nc ^ i 

We next consider the Hubbard chain with the odd 
number of interacting sites. In Fig. 6, the flow of the 
low- lying eigenvalues of Hm/D for Nc = 3 is plotted as 
a function of odd A for several values of the Coulomb in- 
teraction U/{2TTt) = 0.0, 0.5, 1.0, and 1.5. The flow of the 
eigenvalues of Hm /D is quite different from that for even 
Ac- In the case of odd Ac, a number of noninteracting 
sites in the leads are required to reach the flxed point 
that determines the low-energy properties. The number 
of NRG iterations A* that is needed to get the conver- 
gent results increases with U, and in the case of Fig. 6 
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Table II. Low-lying fixed-point eigenvalues of if at for AThub = 3. 



it is estimated to be (b) N* ~ 15, (c) N* ~ 30, and 

(d) N* ~ 50. It means that there is a characteristic en- 
ergy scale determined by T* ~ DA~^^ -i)/2^ where the 
factor A~(^ jg introduced to recover the original 

energy scale from Hn defined in cq. (20). This charac- 
teristic energy scale is determined by width of the Kondo 
resonance, T* Tk, appearing at the Fermi level. 

For > TV*, the fixed-point eigenvalues do not de- 
pend on U and agree with those of the the noninter- 
acting leads determined by eqs. (31) and (32). Namely, 
these low-lying many-body states, E* , have the one-to- 
one correspondence with the quasi-particle states de- 
scribed by eq. (33). The precise correspondence is sum- 
marized in Table II, and these assignments coincide with 
those for the single Anderson impurity. The first and 
third excited states, and E^, correspond to the first 
two one-particle states, the energy of which are given 
by e* = 0.53124 and = 1.42546. We have also con- 
firmed for Nc = 3 that the fixed-point eigenvalues at 
A'' > N* do not depend on whether N is even or odd. 
These feature show that the low-temperature properties 
at T < T* can described by the quasi-particles of the 
local Fermi theory, and it justifies the assumptions made 
in deducing the unitary-limit conductance, gj^^ = 2e^/h, 
at T = for odd Nc- 

To capture the low-energy Kondo behavior at T < Tk 
correctly, one needs to repeat the NRG iterations up to 
N > N* as mentioned in the above. In other words, a 
sufficiently large number of the noninteracting sites are 
required for the reservoirs to make the finite-size energy 
separation smaller than T*, and a similar notice has re- 
cently been emphasized by several authors. '^^'■^^ From 
the results obtained for a small cluster with size cor- 
responding to N < N*, it is still possible to deduce the 
high-temperature properties at T > Tk-^^'^^ In the NRG 
method the logarithmic discretization of the conduction 
band, eq. (20), yields the hopping matrix clement that 
decreases exponentially with increasing TV. 21-23 ^j^^ 
makes the convergence to the fixed point efficient. 

5. Discussion and Summary 

In the present work, we have clarified the difference 
in the transport properties between the Hubbard chain 
of even and odd Nq- Then, what happens in the limit 
of large iVc? It might sound somewhat puzzling since 
the ground state of the Hubbard chain is an insulating 
state in the thermodynamic limit. However, the existence 
of the energy scale T* ~ Tk for odd Nc brings us the 
answer. For odd Nc (= 2M + 1), there is the Kondo 
resonance at the Fermi level. The width Tk decreases 
when Nc increases, and finally Tk ^ in the limit of 



Nc 00 because the Hubbard gap Agap evolving with 
increasing Nc disturbs the electrons to screen the local 
moment. Thus, the value of the conductance g,^ depends 
on the order of taking the limits of Nc 00 and T — > 0. 
For finite Nc, the Kondo behavior at T < Tk causes the 
unitary-limit behavior 



lim 



,].™ 52M+1 



2e^/h . 



(38) 



In contrast, when the thermodynamic limit Nc — > 00 is 
taken first at small but finite T, the conductance must be 
determined by the Mott-Hubbard behavior at T > Tk, 
as 



lim 



lim 

M—KX 



92M+1 







(39) 



Therefore, in order to observe the unitary-limit behavior 
at an accessible temperature, the number of the interact- 
ing sites Nc should not be so large. For even Nc (= 2M) 
there is no Kondo resonance at Fermi level, and the con- 
ductance does not depend on the limiting procedure 



Jim 92M 







(40) 



The two limits, Nc 00 and T ^ 0, considered here 
are analogous to the k ^ and u) ^ limits of the 
vertex corrections for the interacting Fermi systems with 
the translational invariance.'^^ 

In summary, we have studied the conductance through 
a finite Hubbard chain of the size Nc connected to two 
noninteracting leads using the NRG method. The results 
show that the low-lying energy states can be described 
by the quasi-particles of a local Fermi liquid. We have 
also presented a formulation for deducing the dc conduc- 
tance from the fixed-point Hamiltonian. The results of 
the conductance for even Nc show the expected expo- 
nential decay as a function of U at half-filling. 
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Appendix: Reduced matrix element 

In the electron- hole symmetric case, the Wigner- 
Eckart theorem for the spin and axial charge SU (2)spin x 
'S't/'(2)axiai yields 

{I, Iz, Sz ; ^|/n,!ycrl-^ ) -^^) 'S' , 5^ ; r ) 

= {S', S',; 1/2, a\S, S,) {!', I',; 1/2, 1/2|7, /,) 



X Fr,,,{I,S;r\I',S';r') 



(A-l) 



where the Clebsh-Gordan coefficient appears for the 
total spin {S' , S'z;l/2,a\S, Sz) and for the total axial 
charge (/', 1/2, 1/2|/, J^). The invariant matrix ele- 
ment Fn^„{a'\a) has the following properties against the 
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9 



exchange of the arguments a' and a; 
FnAI,S - 1/2- r'\I- 1/2, S;r) 



= (-1)" 



25+1 



21 



2S V 27 + 1 

xFr,,,{I-l/2,S;r\I,S-l/2;r'), 

FnAI,S+l/2;r'\I-l/2,S;r) 



(A-2) 



n+l 



2S+1 / 21 



^ V 25 + 2 V 2/ + 1 

xF„,,(/-l/2,5;r|/,5+l/2;r'). 



(A-3) 
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